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Abstract 

We develop the general theory of stars in Saa's model of gravity with 
propagating torsion and study the basic stationary state of neutron star. 
Our numerical results show that the torsion force decreases the role of 
the gravity in the star configuration leading to significant changes in the 
neutron star masses depending on the equation of state of star matter. The 
inconsistency of the Saa's model with Roll-Krotkov-Dicke and Braginsky- 
Panov experiments is discussed. 
PACS number(s): 04.40.Dg,04.40.-b,04.50.+h 
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1 Introduction 



In recent years the interest in scalar-tensor theories of gravity has been renewed. 
One reason for this is the important role which these theories play in the under- 
standing of inflantionary epoch. On the other hand the scalar-tensor gravitation 
(the so called " dilaton gravity" ) arises naturally from the low-energy limit of the 
super-string theory [0, @. 

The predictions of scalar-tensor theories may differ drastically from these of 
general relativity. For example such a phenomenon - "spontaneous scalarization" 
was recently discovered by Damour and Esposito-Farese as a non-perturbative 
strong field effect in a massive neutron star ||. Other interesting phenomenon 
is "gravitational memory" of black holes proposed in @]. The "gravitational 
memory" in the case of boson stars was investigated in [H] (see also ||). Their 
stability through cosmic history using catastrophe theory was investigated in [[5] . 

Many theories of gravity with propagating torsion involving a scalar field 
have been proposed in the last decades, too ||, [f|, |fL0| . In such theories con- 



trary to the usual Einstein-Cartan gravity ||TTJ-||T3|], there are long-range torsion 



mediated interactions. Carrol and Field [14] have examined some observational 
consequences of propagating torsion in a wide class of models involving a scalar 
field. They conclude that for reasonable models the torsion could be detected 
experimentally. 

Recently a new interesting model with propagating torsion was proposed by 



Saa |15|-||19||. This model involves a non-minimally coupled scalar field as a 
potential of the torsion of space-time. As one can see Saa's model is very close 
to the dilaton gravity. 

In the present article we investigate both analytically and numerically a neu- 
tron star in the Saa's model and compare obtained results with these in the general 
relativity. We also discuss new predictions of the theory under consideration. 

The paper is organized as follows. In section 2 we consider briefly Saa's 
model. In section 3 we give the necessary information for the vacuum solutions 
of the field equations. The equations determining static equilibrium solutions 
for a neutron star are discussed in section 4. Numerical results for the neutron 
star are discussed in section 5. The stability of the neutron star is discussed 
via catastrophe theory in section 6. The inconsistency of the Saa's model with 
Roll-Krotkov-Dicke and Braginsky-Panov experiments is discussed in section 7. 



2 The model with torsion-dilaton field 



Here we give a brief description of Saa's model. For more details one can see 



1§-[|T71, [flj. 

Consider four-dimensional Einstein-Cartan manifold .A/f^ 1 ' 3 **, i.e. four- 
dimensional manifold equipped with metric g a p and affine connection T^ 7 
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with torsion tensor 5 Qj g 7 . 

The main idea of articles [|T^]-[[r7| is to make the volume form d A Vol 
compatible with the affine connection on the Einstein-Cartan manifold Jvi^ 1 ^ 
via the compatibility condition: 

£ v (d A Voi) = iy^)d A Voi, (i) 

where £ v is the Lie derivative along an arbitrary vector field v and V M is the 
covariant derivative with respect to the affine connection. It turns out that com- 
patibility condition (|l]) is fulfilled if and only if the torsion vector 

9 - 1 <? v 

is potential, i.e. if there exists a potential B, such that 

s a = v Q e = d a e. (2) 

In this case Saa's condition fll]) implies the form 

d A Vol = f(x)d A x = e- m ^\d i x (3) 

of the volume element in Einstein-Cartan manifold. As it was pointed out in 
p0| compatibility condition ([]]) leads to covariantly constant scalar density / = 



e~ 36 y \g\ with respect to the transposed connection (T T ) a/ i t = T^ 7 , not with 
respect to the usual connection T^ 7 . Therefore the Einstein-Cartan manifold for 
which compatibility condition (fj) is fulfilled was called transposed-equi-affine and 
the corresponding theory of gravity - transposed-equi-affine theory of gravity. 

The most important mathematical consequence of the condition (|J) which 
leads to new equations of gravity is the generalized Gauss' formula: 

/ d A Vol(V^)= I rf 3 S/. (4) 

J M JdM 

The natural choice of the lagrangian density for gravity is: 

C g = - f K R = - t (r +6V„S" + 12^^ - K^K^A , (5) 

c being the velocity of light, k = 8ttc~ 2 G being the Einstein constant, G being the 
Newton constant. Here R = g a/3 R a p is the scalar curvature with respect to the 
affine connection, K^x = K^ u \ + 2g^ u S\] is the traceless part of the contorsion: 

K^^ IU = K^v^ = 0, and R is the scalar curvature with respect to the Levi-Civita 
connection. 

The traceless part of the torsion doesn't vanish only if spin-non-zero matter 
presents. In the present article we consider only spinless matter (as we know from 
Einstein-Cartan theory of gravity, the effects due to the spin become essential 
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at density over 10 57 g/cm 3 []TT] which is too far from the physics in the stars). 
Therefore we put K a ^ = and obtain a semi-symmetric affine connection: 



Sa/3 1 = S [a 5jy (6) 



In this case we have: 



£ G = -fR = -f K {R +6V„S" + 12S^ . (7) 

Denoting the lagrangian density for the matter by Cm and using the volume 
element (PD we write down the action of gravity and matter in the form: 



g + A m = -%J d A Vol R + i J d A VolC M . 



Due to the new Gauss' formula (||) the term SV^S^ in the lagrangian ([?]) gives a 
surface term in the action integral and doesn't contribute to the equations of 
motion. Hence, these equations may be derived from the modified action: 



(9) 



This action is very close to the one of the dilatonic gravity arising from low- 
energy limit of the superstring theory. Two essential differences between our case 
and the dilatonic one are: 1) the matter action includes the dilaton-like term e~ 3e 
which arises in a natural way, as a part of the volume element of space-time, and 
2) the sign before the term 12S^S^. Following the above described reasons we 
call the field 0, which originates from the space-time torsion and plays the role 
of the dilaton field in Saa's model, "a torsion-dilaton field". 

Taking variations with respect to the metric g a p and torsion-dilaton field G, 
and using the generalized Gauss' formula, we obtain the following equations of 
motion for the geometrical fields g and 9: 



c?^ + v Al v„e-^ne= %T t 



X 8C 



□G = f \C M - l-flfj - \R- (10) 

Here G^ v = R^ u — \gfj, u is the Einstein tensor for the affine connection, its trace is 
G = g^ u G^ v = —R; T^ v = bLul^g^ is the symmetric energy-momentum tensor 
of the matter ; its trace is T = g^T^ and V a S a = ^V^G = DQ. From the 
first equation of the system (|Tt]) it follows that: 

R = -3aQ -^T (11) 



4 



Then combining this result with the second equation of the system ( |T0D we obtain: 

v a s° = ne = -%(c u - \ 5 4^ + It) . (12) 



The equation ( 12|) shows that under proper boundary conditions, and in the pres- 
ence only of spinless matter, the torsion-dilaton field G is completely determined 
by the matter distribution. Further on, as a basic system we will use the system: 



G^ + v„v I/ e-^ne= fr^, 
v a s° = ne = -%(c u - \ 5 -^- + It) . (13) 



From this system one can derive (using Bianchi identity) the differential conse- 
quence: 

V TZ + TZS a = ^RS a (14) 

which is a generalization of the well-known conservation law = in general 

relativity. 

To have a complete set of dynamical equations one has to add to the above 
relations the equations of motion of the very matter. For the purpose of the 
present article we need to consider only a perfect fluid. Its theory was recently 
described in [p0|| . Here we give the basic results. 

The continuity condition describing the conservation of the fluid matter can 
be written in the form: 



L 



d s E a n(x)u a (x) = 0, (15) 

where u a {x) is the fluid four- velocity, normalized by the relation g a pu a u^ = 1, 
n(x) is properly defined a fluid density, d 3 E a is a proper three dimensional surface 
element depending on the choice of the volume element via the Gauss' formula, 
and A^ 1 ' 3 ) is an arbitrary domain. 

Considering the volume element @ as an universal one we must use it in the 
continuity condition, too. Therefore according to the generalized Gauss' formula 
we can rewrite relation (|15D in the form of a continuity equation of autoparallel 
type: 

V Q (n(x)u a (x)) = 0. (16) 

We take the lagrangian of the fluid with internal pressure p in the usual form: 

C ll = —e = —nc 2 — nil, (17) 

where II is the elastic potential energy of the fluid ; (ffl = —pd(^) and the 
symbol " d " means that the corresponding differential form isn't exact. Taking 
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into account the relation — = p it's not difficult to obtain the equations 

of motion for geometrical fields g a p and G in presence of perfect fluid: 



C?^ + (V M V„-^n)e = ^(£ + P )«Y-pir 

□9 = -£(e-p); (18) 

In addition one can show that: 

V CT T Q CT = {e + p)(8 a a -u a u a )S <n or 

V«T° = 3{e+p)u°u a S a . (19) 

Making use of (|T9|) and of the continuity condition (]T6|) one can obtain the equa- 
tions of motion of the perfect fluid (just as in general relativity): 

(e + p)u p V p u a = (5* - u a u?) VpP- (20) 

The equations (|19f) are equations of a geodesic type. In particular, considering 
dust matter (p = 0) we have: 

u p V P u a = 0, (21) 

i.e. we can conclude (just as in general relativity) that a test particle in the theory 
under consideration will move on a geodesic line. We will need this conclusion in 
the next sections. For more details concerning the relativistic perfect fluid in the 



theory under consideration we refer to [20]. 



3 Spherically symmetric vacuum solution 

The asymptotic flat, static and spherically symmetric general solutions of the 
vacuum geometrical field equations (|T0|) are known [pi], f22fl . In Schwarzschild's 
coordinates they are described as a two parameter - {K, a} family of solutions 1 : 



9oo = e , 

r = ^ae 2 sinh ( 2 ) > 
We use asymptotic conditions v ^ 0, 0-^Oatr-^oo without loss of generality. 
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3 (if - |) 2 + I, 



S r = 0' = K£. (22) 

Here and further on the prime denotes a differentiation with respect to the vari- 
able r. All quantities in formulae (p2|) are represented as functions of the variable 
v. This is the most convenient form of the vacuum solutions. 

The parameter K presents the ratio of the torsion force (as defined in f20|| ) 
and the gravitational one: K = S r /( \v'\ In the case when K = we have the 
usual torsionless Schwarzschild's solution and a = r g is the standard gravitational 
radius r g . 

In the model under consideration the value of the fundamental parameter of 
the theory K (which is constant in vacuum) is not an independent integration 
constant. Instead, we shall show that it is determined by the total mass of the 
star, or by its radius and depends on the matter distribution, on the equation of 
state of the star's matter, and so on via the solution of the full system of equation 
of the star's state. 

The parameter a is positive (a > 0), and may take arbitrary values. It is 
related to the total mass of the star, too. 

The asymptotic behaviour of the solution is: 

2M K , G 

goo „ i 23 

err 

2G(M Kepler - M e ) 

9u ~ 1 + 1 , (24) 

c r 

where M Kepler = \%- p describes the asymptotic dependence of goo on the variable 

r, and the mass Mq = 3|r lim^oo r 2 S r describes the asymptotic dependence of 
the torsion-dilaton potential on r (or the asymptotic of S r ). The mass Me may 
be considered as a source of torsion-dilaton force and is analogous to the "scalar 
mass" introduced in [ 23| . As one can see the scalar mass Mq depends on K. 



In the model under consideration, a test particle moves along geodesic lines. 
Therefore, the keplerian mass measured by a test particle in the asymptotic region 
of space-time will be the mass M Kepler = ^ lim r ^oo \r 2 v'e v . As we will see in the 
next section, the mass M Kep[er is positively defined. The appearance of two masses 
in the asymptotic solution is related to the violation of the strong equivalence 
principle (the weak equivalence principle for test particles is not violated). In our 
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case the ratio M e to M Kepler is just 3K and depends on the ratio K of the torsion- 
dilaton force to gravitational one in vacuum. As in Brans-Dicke theory one can 
define so called tensor mass Mt, which is the mass measured by a test particle 
in Einstein frame, i.e. a test particle moving along a geodesic in space-time with 
metric g^ v = e~ 3@ g^ u . It's not difficult to see that Mt and M Kepler are related by 
the formula: 

M T =(l-lK]M Kevler . (25) 
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In the next section we show that the tensor mass My is also positively defined. 
The existence of two positively defined masses M Kep i er and M T makes the question 
about the energy of the geometric-field complex g, subtle - which of them is 
the true measure of the energy of the star? 

It turns out that one has to choose the tensor mass Mt- The reasons for this 
we will comment shortly in the next section. For detailed consideration of the 
question of the star mass definition in presence of a scalar field we refer the reader 



to 24 . 



4 The Basic equations for a star 
4.1 General considerations 

Here we will discuss some general properties of the system of equations of the 
star without specifying the matter's equation of state. 



The system of the equation ( |18| ) for geometrical fields g and can be rewrit- 
ten in the form: 

§ <T S ff -3S <T S* = -%(e-p), 

I£ = (e + pKu M -pg£, (26) 

where V is the covariant derivative with respect to the Levi-Civita connection, 
G 1 ^ is the corresponding Einstein's tensor and 

^ = V»S V + S^S V - g'fivS* + fiSrS*. (27) 

In this paper we restrict ourselves with consideration of the static and spher- 
ically symmetric case. Hence, the metric has the form 

ds 2 = e v(r \cdtf - e x{r) dr 2 - r 2 (d6 2 + sin 2 (#) dip 2 ) (28) 
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in which the functions v = u(r), A = A(r) depend only on the Schwarzschild's 
radial coordinate r and the torsion-dilaton field depends only on r, too. 

In this case we obtain the following equations for the functions u, A, G, and p: 



A ( r 2 r ) r 2 ~~ c? £ ^Sq, 

( ~ + ~?) + ^ = — ^P ~~ 



-e 



e- A (# + ( + 2 ) S r - 3S r 2 ) = £(e - p), 

p' = - \(e + p)u\ 

p = p(e). (29) 

Here p = p(e) is the matter's equation of state and correspondingly: 



— ( \v S r + ^S r — 2S r J e , 

£2 = ^/ + l u > Sr _ l X t Sr + 1^ _ ^ e -A_ (30) 

{} 

The second equation of the system ( p9]) (i.e. G 1 ! = ^T/ ) may be consid- 

ered as a constraint creating a relation between SV, v' and e A , namely: 

, 1 + ru' -QrS r -lr 2 u'S r + 6r 2 S r 2 , . 

* x TTJ^ ' (31) 

Using this relation we can put our system ( p9|) in a normal form: 



Q' = S r , 



5;. = - ^ + (f(£-p)- f(e-p)rS r )e A , 

p' = -(e+p)f, 
p = p(e), 



A _ 1 + 2r£ - 6rS r 

1+ fpr 2 

It's seen that the first two equations are separated, and the rest generate a sub- 
system independent of them. 

The equations fl32|) must be solved with proper initial and boundary condi- 
tions. From a physical point of view the solutions regular at the center are the 



3r 2 £S r + 6r 2 S r 
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most interesting. The regularity means that there exists a local lorentzian system 
in neighborhood of the center, i.e. e A (°) = 1 and the pressure is finite at r = 0. 
Hence, we have lim r ^ r£(r) = 0, otherwise, as it may be seen from equation for 
p, the pressure will have at least logarithmic singularity at the center. On the 
other hand the condition e A( ^ = 1 requires that lim r ^ rS r (r) = 0, too. The 
expansion of the equations for £ and S r around the center is: 

r = 

r 

S' r = (33) 

Hence, the behaviour of £ and S r around r = is cons \ ant . To fulfill the above 
restrictions at r — > we must put constant = 0. Hence, we obtain £(0) = 0, 
S r (0) =0. As a final result these considerations imply the following initial 
conditions: 

£(0)=0, S r (0) = 0, P (0)=Pc, (e(0)=e c ). (34) 

At the star surface r = R we have to match interior solution with the exterior 
(vacuum) solution. We will consider the model of the star without surface tension, 
hence p(R) = 0. Then the matter distribution must be continuous at the surface 
of the star and one can show that £ and S r must be continuous at r = R. 
Obviously v and O must be continuous at the star surface, too. Using matching 
conditions: 

S r (R) = S r ext (R), (35) 

we can obtain the vacuum solutions parameters K and a as functions of e c , i.e. 

K = K(e c ) ) a = a(e c ). (36) 

For arbitrary values v c = z/(0) and O c = O(0), v and O will not fulfill the 
matching conditions: 

V [R) = u ext {R), 

Q(R) = Q ext (R). (37) 

Therefore, the separated equations v' — 2£ and O' = S r must be solved under 
proper initial conditions in the following form: 



rR 

v c = v ext 

10 



(R) - f 2£(r) dr, 
Jo 

© c = Q ext (R) - f e( r )dr. (38) 
Jo 
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As a result we obtain all parameters K, a, v cy C , R as functions only of the 
central density e c . Hence, the whole geometry of the space-time in vacuum and 
in the star is completely determined by the matter which carries only the same 
properties described by mass, matter density, pressure, equation of state and so 
on which are familiar from the general relativity. A very important feature of 
the model under consideration is that we are not forced to assign to the matter 
new properties, charges, or something else. Nevertheless we have a new geometric 
field (the torsion-dilaton field G) the physical problem is well defined by the usual 
physical properties of the matter. 

Let's go back to the subsystem: 

£' = -£+(fe- J" £(e-pK)e\ 
S' r = - f + ($(e-p)- %(e-p)rSr)e x , 

p' = -( £ +p)€> 

P = p{e), 



A 1 + 2r£ - 6rS r - 3r 2 £S r + 6r 2 S r 



e A = ^ 7 ZJ ■ (39) 



We can't define a local gravitational mass in the form mGn{r) = 7^r(l— e~ A ) , 
as in general relativity because in our case mGn{r) is in general not positively 
defined (see for example the vacuum solution). In the theory under consideration 
we define the local mass as m v (r) = ^r 2 £(r) . This local mass is connected with 
local Keplerian mass Tn Kepleir (r) (^ ne mass measured by a non-self gravitating test 



particle in a circular geodesic orbit with radius r [24]) by the relation 

m K, v iS r ) = m u (r)e u . 
Similarly, we can define a local scalar mass mj(r) = 3^r 2 S r (r) . We can intro- 



duce also a local tensor mass mr{r) V% as 



(m u (r) - |m e J e l 

m T (r) = ± ? — — 

(1 _ Gm & ) 



The initial conditions imply m u ( y 0) = m K kr (0) = mg(0) = mj'(O) = 0. 

Now the system ( p9[) can be rewritten in terms of the masses m v {f) and me(r): 

I ( -\ Mi ft I r~ ^,2 ^ „\\ m„ | I6ir^2^„\ 

v 



m' u = (!-(!+ £(e - p)r A ) e A j ^ + ^r 2 se 
m' e =(l-(l+ ^(s-p)r 2 )e x ) ^ + ^ r 2 {s-p)e\ 

p' = -^(e + P ) 

P = p{e), 
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This system is a generalization of Tolman-Oppenheimer-Volkoff's one for a star 
in general relativity p5| . Using the first and the second equation, it's not difficult 
to show that m u (r) and Tng(r) are positively defined. Indeed, taking into account 
regularity at the center we obtain: 



e A(r) f r e -A(r) Iftr^ 

Jo c 

me = e A{r) f e~ A{r) ^(e - p)r 2 dr, (41) 
Jo c 



rr 



where 



A(r) = r ^^ 2 y d r. 

Jo r 

In the same way from the above equations we have: 

K - m e )' = (l - (l + £ (e - p)r 2 ) e x ) - ^r 2 (e - 3p)e\ (42) 

Solving this equation with an initial condition [m v — mg)(0) = we obtain : 



r 



m u - me = - e Mr) / e -A(r)/ £ _ 3p ) dr _ ( 43 ) 

JO 

Hence, it's seen that m v — < inside and outside the matter if e — 3p > 0. 
In other words we obtain for k(r) = \ m9 ^\ that k > 4 and k takes a value \ when 

v / 3 m v (r) — 3 3 

the matter is ultrarelativistic (e = 3p). The parameter k takes its maximum 
value | in the case of nonrelativistic matter (e ^> p). If we assume following 
Zel'dovich p9 |, [p0| that e < 3p may happen, then in general the vacuum value 
of k which is just K = k(R) may change its sign passing through the zero at 
e = p. For realistic equations of state we obtain K 6 [|, |]. Therefore we have 
M T e[lM Kepler , \M KEpl J (see©). 

For completeness we will give an expression which is a generalization of the 
well-known Tolman's formula [32]. From equations (|26| ) we have 

i?° + 3Vo5°= f (T °- IT) - fne. (44) 

In the static and spherically symmetric case, one may show that the following 
relation is fulfilled: 

R°o + 3 Vo S° = —^i=d a ( vfe- 3 V4Vl , (45) 
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where T a (P = {^} are Christoffel symbols. Hence, we obtain 



/ (r° + 3 Vo S°) e~ 3 °^\d 3 x = I g^ { U P a e-™^\dll a = 4tt %^.(46) 



Therefore for Keplerian mass we can write down 



,.2 



M Kepler = & J M \Ro + 3 Vo S°\ e-^^x 



$j (2T°-T)e- 3e ^\d 3 x- £e I J UQe^^x. (47) 



Taking into account that 

-30 /i I _ 4tt M e G 



/ □ee- 30 ^jd 3 a: 



3 c 2 



and 2Tq — T = e + 3p we obtain: 



A epler 



I M 

On the other hand, taking into account that 



h [ (e + 3p)e- 3e J\^\d 3 x+ \M e . (48) 

JM 



Me = 3KM Kepler = % f (e - p)e~ 3e ^\d 

JM 



one can rewrite the above formula in the form 



' 3 x 



M K , 

K epler 



^ i / (e + 3p)e- 3e ^\d 3 x 



f ee~ 3e Md 3 x= ±£ [ R ee^-^dr. (49) 



.M "JO 



Comparing this formula with (^) we immediately obtain the following expression 
for the tensor mass 

M T = \ [ (s + 3p) e - 3e v& 3 z = % \* (e + 3p)e^- 3e r 2 dr. (50) 

CT JM V C JO 

From the explicit expressions for the masses we have 

< M T < M Kepler < M (51) 
when matter satisfies the condition e — 3p > 0. 

The number of the particles in the theory under consideration is given by 

N = 4vr / n(r)e~ 30 r 2 dr (52) 
Jo 
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where n(r) is the particle density. We note that we consider cold matter and we 
have e = e(n) , p = p(n). Therefore for the rest mass Mr we have 



R 



M R = m N N = 4irm N / n(r)e 2 r dr (53) 

Jo 

where is the nucleon mass. 

The binding energy of the star Eb depends on the choice of the star mass. 



As it is shown in |24| the right choice must satisfies the following physically 
reasonable conditions in the spherically symmetric case: 

1) to be non-negative; 

2) to be identically zero in Minkowski space-time; 

3) to be non-decreasing function of r; 

4) the maximum of the mass to coincide with the maximum of the rest mass( 
the number of the particles). 

As it is seen from numerical calculations both Keplerian and tensor mass 
satisfy the first three conditions but only the tensor mass satisfies the fourth. 
Therefore binding energy must be calculated with respect to the tensor mass: 

Eb = M T - M R . (54) 
4.2 Neutron star model 



First we consider non- interacting neutron gas at zero temperature | 28| , [|33 1 . The 



energy density, the pressure and the particle number density in a proper normal- 
ization are given by: 



4 5 

m N c 



3vr 2 /i 3 



4 5 

m N c 3 , , 

71 = 3^ (55) 



where 



g{v) = s S/V/^ + Z^ ~ V^ + /i 2 (2/i-3) -3in(v^+ Jl + fi) 



M = + ^ 2 (2/i - 3) + 3/72(7/1 +y /l + f i)), (56) 

^ = ( ^^Tlf ) 2 ' iFermi being the Fermi's momentum, mjy being the neutron mass. 

We are interested in the difference between the predictions of the theory under 
consideration and of the general relativity. For this purpose the equation of state 
for a non- interacting neutron gas is sufficient. 
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As a more realistic equation of state we consider the analytical approximation 
(according to Zel'dovich and Novikov P0| ) of Tsuruta-Cameron's equation of 
state ||31|| . In this case the interaction between the nucleons is taken into account 
in a simple approximation and the pressure is given by: 

p c 2 p c 2 f As \ 1/2 
P = £ + — - — { 1 + ^) (57) 

where po = 5 * 10 15 g/cm 3 as the relation between the particle number density 
and the energy density is 

e = itlnc 2 I n + tun — ] (58) 
V PoJ 

It turns out that these two examples present typical results which qualitatively 
agree with the results for the other equations of state of star's matter. 
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5 Numerical results and discussions 



We have solved the system of equations (^) coupled with the state equations 
( |56| ) and fl5T|) numerically using the method due to Runge-Kutta-Merson with 
automatic error control. The results are shown in the corresponding figures. 
Hereafter all masses are measured in units Mq. 
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Figure 1: a) M — log(e c ) dependence. 
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b) Mt — Mr dependence. 



First we concentrate our attention on the case of non-interacting neutron gas. 
In Fig. la) the dependence of the three masses Mt,Mk,Mr on the central density 
e c is shown. The appearance of a cusp in Fig. lb) , where the dependence of the 
tensor mass Mj- on the rest mass Mr is presented, shows that their maxima lie 
at the same point. Although, the maxima of the rest mass Mr and Keplerian 
mass Mk are too close in Fig. la) they don't lie at the same point, as it may be 
seen from Fig. 2a) which shows the dependence of Mk on Mr. We see also that 
Keplerian mass is considerably greater than the tensor one - about three times. 

In Fig. 3a) the My — R dependence is represented. It's seen that the Mt — R 
curve in our case is fairly similar to the one of general relativity, but there are 
significant differences, too. The maximum mass Mxmax i n our case is ~ 0.35Mq, 
while in general relativity the Oppenheimer-Volkoff's mass is Mov — 0.7Mq. 
The radius corresponding to the mass Mq is R = 4.2km, while in the case of 
general relativity R = 9.6km. If we look at Fig. 2b) where the dependence of Mt 
on the central density e c is shown, we note that MT max lies at e c ~ 4.5*10 16 g/cm 3 , 
while Mov lies & t £ c ~ 5 * 10 15 g/cm 3 in general relativity. The average density 
in our case is about 4 times greater than the one in general relativity. Hence, in 
the model under consideration the neutron star is more compact and has a mass 
about 1/2 - times smaller than Mov- In Fig. 3b) the dependence of Keplerian 
mass on the star radius is presented. It's seen that the Keplerian mass is about 
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1.5 times greater than M ov . 
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Figure 2: a) M Kepler — M R dependence. 
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b) M T — log(e c ) dependence. 
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Figure 3: a) M T — R dependence. 
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b) Mxepier — R dependence. 



In the Fig. 4a) the dependence k(r) is shown inside the star (for central density 
7.5* 10 15 (?/cm 3 ). In accordance to the general considerations k increases from the 
center of the star to the surface, where k takes a value K = k(R) ss 0.45 — 0.46, 
which is close to 0.5. The dependence K(e c ) of K on the star central density e c 
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Figure 4: a) k — r dependence. 
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b) K — log(e c ) dependence. 



is shown in Fig. 4b). It's seen that K decreases when density increases, which is 
similar to the previous case. So, the ratio of the torsion force to the gravitational 
one takes its minimum value at the center of the star and is the greatest at the 
surface. 




As it may be seen from Fig. 5a) expressing the dependence K(Mt), the 
torsion-urged effects are relatively strongest in the case of small masses - with 
increasing of the star mass (up to the point where the star loses its stability) K 
decreases. It's seen from Fig. 5b), where the dependence of K on the star radius 
R is shown, that in the area of stability K decreases when R decreases too - the 
more compact stars are, the smaller K they have. 

Fig. 6a) presents the dependencies 6{r) and v{r) inside the star. One may 
see that \v — 30 < everywhere. The dependencies m T (r) ,m Kep i er (r) and mg(r) 
are shown in Fig. 6b) for central density 7.5 * 10 15 . As it has already mentioned 
all masses increase with r. 
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b)rriT, m Kepler i m e — r dependence. 



The following figures illustrate the case of Tsuruta-Cameron equation of state 
(TCES). 
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Figure 7: TCES. a) M - log(e c ) dependence . 



b) M T — M R dependence. 



We see from the figures that the maximum tensor mass in this case is about 
2Mq and the corresponding radius is about 7.5km - the same quantities in 
general relativity are correspondingly ps 1.6Mq and m 11.5km. Hence, the 
interaction between the nucleons leads to an increase in the maximum mass, as 
in general relativity. 

Note the differences between the Fig. 5b) and Fig. 6b) (for the case of non- 
interacting neutron gas), and the corresponding Fig. lib) and Fig. 12b) (for 
the case of Tsuruta-Cameron equation of state). There one can see the strong 
dependence of some results in the Saa's model of gravity with propagating torsion 
on the equation of state of star's matter. 
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Figure 8: TCES. a) M Kepler — Mr dependence. b) Mr — log(e c ) dependence. 
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Figure 9: TCES. a) M T — R dependence. 



b) M Kep i er — R dependence. 



We have also examined the Harrison- Wheeler's equation of state p5 |. As in 



general relativity the numerical results are very close to these for the noninter- 
acting neutron gas. For example the maximum tensor and Keplerian mass is 
correspondingly « 0.35Mq and ~ IMq, and the corresponding radius is 3.8km. 

Other equations of state (of politropic type) have been examined, too. The 
corresponding maximum tensor mass of a neutron star reaches a value about 
2.5 — 2.6Mq, while the corresponding Keplerian mass is about 6 — 6.5Mq. 

As it is seen from numerical calculations the tensor mass and Keplerian mass 
differ very significantly from each other. This behaviour of the masses is qualita- 
tively the same as in the case of a boson star in Brans-Dicke theory with u — — 1 
|24|. Saa's model corresponds to the value of the Brans-Dicke parameter uj = — | 
which is close to —1. That's why the observed qualitative agreement is natural. 
Hence, including a scalar field with an approximately the same uj in physically 
different kinds of stars we find similar departures of the corresponding predictions 
of general relativity. 

This conclusion holds also in the case of large uj. For example, in || and 
Brans-Dicke boson star has been considered with uj = 400. There the results are 
quite similar to corresponding ones in general relativity which is just the limit 
uj — > oo. The star is a bit lighter but has a higher density in future. 
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Figure 10: TCES. a) k - r dependence. 



b) K — log(e c ) dependence. 




Figure 11: TCES. a) K - M T dependence. 
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6 Stability analysis via catastrophe theory 



Methods of catastrophe theory have recently been applied for investigation the 



stability of self-gravitating systems as neutron and boson stars in [2!| and 
(in the case of boson stars see also [§] and [[j]). 

Here we discuss briefly the important question of the stability of a neutron 
star in our case using tools of catastrophe theory 

The basis for the stability analysis are figures Fig. lb and Fig. 7b. The con- 
served quantities which we use to calculate the binding energy are the tensor 
mass and the rest mass (particle number). Drawing them against each other we 
obtain so-called bifurcation diagram. Actually, Fig. lb and Fig. 7b are bifurcation 
diagrams. When the central density increases one meets a cusp. The appearance 
of a cusp itself isn't enough to conclude that the stability of the neutron star 
changes. However, if it's known that the star is stable along the first brunch 
then it will become unstable along the second brunch. In general relativity the 
stability of the neutron stars for small central densities is proved by using linear 
perturbation analysis. As it's seen from the figures, the behaviour of the basic 
conserved quantities in our case is qualitatively the same as in general relativity. 
That's why it's natural to assume that for small central densities, i.e. along the 
first brunch the star in our case is stable against small radial perturbations. Be- 
sides of that the first brunch is composed completely of negative binding energy 
states which shows that it is potentially stable. 

Making use of the above considerations we can conclude that the neutron stars 
with central densities from small values to the cusp are stable. Beyond the cusp 
one radial perturbation mode develops instability and the star becomes unstable. 
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7 Saa's model, Roll-Krotkov-Dicke, Braginsky- 



Panov and Eotvosh experiments 

Despite the obvious beauty of Saa's model we shall show that it contradicts to the 
experiments by Roll-Krotkov-Dicke, Braginsky-Panov and Eotvosh. To do this we 
need the equation of motion of an isolated test body corresponding to the model 
under consideration. In our case the local conservation law for energy-momentum 
tensor can be written in the form (see ( ^4|) and (|TT|)) 



{} 



V CT (e" 3e T a CT 



3e 



-30 



c 



M 




(59) 



Defining the four-momentum and the mass-center of a small test body as 
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and following the same procedure as in general relativity J3"3 |, [ 3"5 ] we obtain the 
equation of motion of such body in our case: 
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We consider a static spherically-symmetric case and therefore we put S° = 0. 
Fixing a local inertial frame in which the test body is at rest at the moment 
considered and the Christoffel symbols vanish at the location of body we have 



d 2 XP 



(61) 



where M is the inertial mass of the small test body. 

Therefore in the model under consideration test bodies with different A would 
accelerate differently in the solar field (or in the earth's field). In good approxi- 
mation we can consider the test body as built by electromagnet ically interacting 
particles in gravitational and torsion field, i.e. the Lagrangian density for the test 
body is taken in the form 



C 



M 



nem , f*pf 
*-M ~r 



(62) 
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where 

rem _ pfiv rp 

~ l 67 r 



and 

(for Lagrangian density for particles in Saa's model see [0). It can be shown 
that 

C M-l-ji = °- (63) 



Consequently we obtain 



A = 3 / £ 
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(64) 



Here £ e and £ m are correspondingly the electric and the magnetic energy containt 
of the test body. Therefore from (pH) we have 



r\ 2 Yi 3 P — P 

-4^r = -3S^X)^—p^. 
dt 2 M 



(65) 



For aluminum and platinum the magnetic energy is small compared to the electric 
one. The ratio ff for aluminum and platinum is |35| 



MJ 



1.7 *1(T 3 , 
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(66) 



£ \ 

= 4.5*10~ 3 . (67) 



On Earth we have 
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Si = 8.4 * KT 3 ^, 



% runs from 1 to 3. 
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Out of the source the gravitational potential and the torsion vector are related 
by Si = KU,i (see (fZ^)) and therefore 



d 2 X,\ d 2 X, 



J Al V ^ 2 J Pt 



8A*10~ 3 KU ti . (69) 



According to the experiments of Roll-Krotkov-Dicke [36] and Braginsky-Panov 



[37j the acceleration of aluminum and platinum don't differ by 1 part of 10 11 or 



10 of U } i respectively in the solar gravitational field. Taking into account that 
| < K < | we conclude that Saa's model contradicts to these experiments. 

It's not difficult to see that Saa's model contradicts to the experimental data in 
Eotvos experiment, too. Indeed, according to Eotvos experiment the acceleration 
of the bodies of different materials in Earth's gravitation field don't differ by 1 
part of 10 11 of g |42| (or 10 12 according to [Q) . Using again relation (^) and 
inequality | < K < | we see that Saa's model contradicts to Eotvos experiment. 



8 Summary 

In this article we have examined the basic spherically symmetric stationary state 
of stars in the Saa's model of gravity with propagating torsion. 

In the model under investigation there is no need to consider unknown charges 
creating the torsion-dilaton field. Its source is the very spinless matter. The whole 
geometry of the space-time (including metric and torsion) is determined by the 
familiar properties of this matter. 

The parameters of the vacuum solution are determined only by the spinless 
matter without adopting an existence of new properties, too. In contrast to the 
corresponding models in the general relativity here we have two parameters K and 
a of the vacuum solutions. The values of these parameters depend on the mass 
distribution in the star which is related with the equation of state of the star's 
matter. For a fixed equation of state both parameters become functions only of 
the star mass, but these functions are not the same for the different equations of 
state. The first parameter K being the ratio of the magnitude of torsion-dilaton 
force and of the magnitude of gravitational force for realistic equations of matter 
state takes values in the interval [|, |], depending on the star's mass. The second 
one - a is analogous to the gravitational radius in general relativity and takes 
positive values depending on the value of the parameter K and on the value of 
the star's mass. 

To be specific in the present article we restrict our attention to the model of 
neutron stars where the effects of nonlinearity are essential as in general relativity. 
Numerical results and analytical considerations show that the space-time torsion 
may have a significant role in their structure. The new torsion force decrease in 
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some extent the role of the gravity in the star configuration and may lead to an 
increasing or decreasing of the maximum neutron star mass depending on the 
equation of state. 

The complete investigation of the consistence of the whole Saa's model of 
gravity with propagating torsion (including all type of physical fields) with the 
reality is still an open problem. The results of the present article may have not 
only independent value, but are necessary for reaching the solution of this critical 
problem. For example, after the first version of the present article was send 
for publication a new results based on it which show the inconsistency of Saa's 
model with solar system gravitational experiments were found and published 
independently fl38| . 

Saa's model is a simple model based on pure geometrical reasons which allows 
to overcome a basic difficulties of the old models of gravity with torsion: the in- 
consistency of the application of the minimal coupling principle in action principle 
and directly in the equations of motion |15 |. Moreover, it leads to propagating 
torsion which is another important physical property of this model. Unfortu- 
nately, it is not compatible with basic physical experiments. This means that one 
has to modify this model preserving its important new physical properties in a 
proper way to comply with real physics. A new investigations in this direction 
are in progress (see for example P5| , PI). 

Nevertheless the model under consideration contradicts to basic experiments 
its detailed investigation is quite instructive because Saa's model is a special 
case of scalar-tensor theory of gravity with a nonminimal coupling of the scalar 
torsion-dilaton field with matter. A similar, but more general coupling of string- 
dilaton can be expected in the string theory and it leads different physical effects. 
For example, in the article by Damour and Polyakov as a consequences of 
nonminimal coupling between matter and dilaton was derived a violation of the 
equivalence principle, as far as some string-dilaton effects in the early universe 
£T|. Unfortunately, the present days string theory is not able to predict definitely 



the form of the coupling between string-dilaton and the real matter. In contrast, 
Saa's model is the only one we know with complete determined interaction of 
the (torsion) dilaton with all kinds of matter. This interaction is similar to the 
one expected in other models. In the present article it is shown at first that the 
nonminimal coupling of the dilaton field with the matter will emerge in extreme 
conditions in a neutron star and that it leads to clear new physical phenomena 
of different kind, some details of which may depend on the equation of matter 
state. Most probably similar effects will appear in other possible modifications 
of the theory of torsion-dilaton, as far as in the other types of theories of dilaton. 
For example the neutron stars may give us a way to a real physics in string 
theories if the string-dilaton interactions with real matter will be established. 
Hence, the most important conclusion of present article is that looking for physical 
manifestation of the dilaton one has to investigate in details its influence on the 
neutron star structure. 
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